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ABSTRACT 

The Hubble Ultra Deep Field source JD2 presented in Mobasher et al. (2005) 
is an unusual galaxy that is very faint at all wavelengths shortward of 1.1/im. 
Photometric redshift fits to data at 0.4 to 8 fim yield a significant probability 
that it is an extremely massive galaxy at z ~ 6.5. In this paper we present new 
photometry at 16 fim and 22 fim from Spitzer Infrared Spectrograph (IRS) peak- 
up imaging of the Great Observatories Origins Deep Survey (GOODS) fields. 
We find that the spectral energy distribution shows a factor of ~4 rise in flux 
density between the 16 /xm and 22 /im bandpass which is most likely due to the 
entrance of polycyclic aromatic hydrocarbon emission features into the 22 /im and 
24/im passbands. The flux ratio between these bandpasses can be best fit by a 
z ~ 1.7 luminous infrared galaxy with a bolometric luminosity of (2— 6) xlO^^ Lq 
corresponding to a star- formation rate of 80 Mq yr^^. The predicted flux den- 
sity values at other longer wavelengths are below the detection limits of current 
instrumentation but such sources could potentially be detected in lensed submil- 
limeter surveys. Re-evaluation of the optical/near-infrared photometry continues 
to favor z > 6 photometric redshift solutions, but we argue that the consistency 
of the multiwavelength parameters of this galaxy with other dusty starbursts fa- 
vor the z ^ 2 mid-infrared photometric redshift. The data presented here provide 
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evidence that optically undetected near-infrared sources which are detected at 
24 yum are most likely dusty, starburst galaxies at a redshift of z ^ 2 with stellar 
masses >10^° Mq. 

Subject headings: cosmology: observations — early universe — galaxies:evolution 
- galaxies:individual (HUDF-JD2) 



1. Introduction 



JD2 is an interes ting object in the Hubble Ultradeep Field (UDF: [Thompson et al.ll2005 



Beckwith et al.ll2006l ). Its non-detection shortward of the 1.1/xm bandpass, even in ultradeep 
Advanced Camera for Surveys (ACS) UDF optical data which are sensitive down to ~29 
AB mag, has been interpreted as being due to the redshifted Lyman break. In addition, it 
displays a break in its broad band spectral energy distribution (SE P) between 2.2um and 
3.6/xm which has been interpreted as the redshifted Balmer break. iMobasher et al.l (120051 ) 
argue that it may be a very massive (6x10^^ Mq) galaxy with very little ongoing star- 
formation at 2 ~ 6.5. The large stellar mass of this galaxy derived from fitting population 
synthesis models to the multiband photometry, and high source density implied by the 
presence of one such object in the small area subtended by the Near-infrared Camera and 
Multi-object Spectrograph (NICMOS) UDF, suggest that UV-faint objects, if indeed at 
such high redshift s, contribute as much to the stellar mass density at 2; ~ 6.5 as rest-frame 
ultraviolet-bright Lyman break galaxies. Since the age of the stellar population in JD2 is 
thought to be ~600 Myr, galaxies like JD2 might harbor remnants of the first epoch of star- 
formation, and play an important role in the reionization of the intergalactic medium (IGM). 
Such an object, if truly at high redshift, would challenge current models of galaxy formation, 
which do not pred ict a large number density of galaxies more massive than 10^^ Mq at z > 6 
JPave et allbooeh . 



The photometric redshift solutions in iMobasher et al.l (120051 ) favor a z ~ 6.5 solution 
but provide a 15% probability that the source is at z < 5, with the most likely alternative 
solution implying a dusty galaxy at 2; ~ 2.5. There are a few inconsistencies in the high- 
redshift interpretation of JD2. The optical/near- infrared SED is best fit by a template 
without significant extinction or star-formation. However, the object is detected in the Great 
Observatories Origins Deep Survey (GOODS) Spitzer Multiband Imaging Photometer and 
Spectrometer (MIPS) 24/im survey with a flux density of 51.4±4/iJy. Since the 24/im flux 
would be dominated by redshifted hot dust and polycyclic aromatic hydrocarbon emission, it 
indicates the presence of dust within t he galaxy. The age -extinction degeneracy is very well 
known for red galaxies. Most recently, IStern et al.l (120061 ) have shown that the optical/near- 
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infrared photometry of dusty, extremely red objects (EROs) is almost identical to that of 
evolved, passive EROs, and that only the detection of reprocessed emission at mid- and 
far-infrared wavelengths can break the degeneracy. Nevertheless, Mobasher et al. (2006) 
suggest that the 24/im emission could be explained by the presence of an obscured active 
nucleus. While this is possible, at 2 ~ 6.5, the Chandra 2—8 keV band traces 15—60 keV 
emission, energies which are relatively immune to absorption. The Chandra detection limit 
corresponds to Lx of 3x10"^^ erg s~^ at 2; = 6.5 while the mid-infrared detection implies vl^y 
at rest 3 /z m of 3x10^^ erg s~^. The 10 /im to X-ray luminosity ratios of local Seyfert nuclei 



are ~3 (e.g. iKrabbe et al.ll200ll ). and the flat spectrum of Seyferts in vLu, implies that either 



large (Nh >> 10^^ cm ^) column densities of neutral gas must be obscuring the hard X-ray 
emission or the active galactic nucleus (AGN) is an unusually low luminosity X-ray source. 



redshift fits are skewed in the Mobasher et al. 



In a recent paper, iDunlop et al.l (120071) qu estion the optical limits in the B, V, i and 
z passbands adopted by iMobasher et al.l (120051) and demonstrate that the values in the 



(I2OO5I ) analysis by adopting optical flux limits 



which are too stringent. By refitting the source with revised optical limits, they conclude 
the source to be at ^ ~ 2.2 rather than ai z ^ 6.5. We also note that JD2 has been visited 



in work by lYan et al.l (12004 ) and lChen fc Marzkel (12004 ) prior to the lMobasher et all (120051 ) 
analysis, and had been classified as a z ~ 3 ERO. 

In this paper, we present new 16 /xm and 22/im photometry of this source from Sp itzer 
Infrared Spectrograph (IRS) peak-up imaging of the GOODS fields (ITeplitz et al.ll2006l ). We 
fit the photometry at these two wavelengths and the published 24 /xm data with a variety of 
dust spectral energy distributions and demonstrate that the mid-infrared photometry is most 
consistent with a redshift of 2; ~ 1.7 for JD2. We assess the accuracy of the optical/near- 
infrared photometry and refit the data with this redshift constraint to determine physical 
properties of the galaxy. We conclude that it is a luminous infrared galaxy with Ljr ~ 
10^^ '' L0 and not a z > 6 galaxy with a dust obscured AGN. Throughout this paper, we 
adopt a Qm = 0.27, = 0.73, Ho=71 km s~^ Mpc~^ cosmology. 



2. IRS Peak-Up Imaging of GOODS-S 



The 16/im and 22/im observations were carried ou t using the peak-u p imaging capa- 
bility of the Infrared Spectrograph (IRS) instrument (IHouck et al.l |2004 ) on the Spitzer 
Space Telescope. The field of view of the peak-up imaging camera is 54"x80". The central 
~130 arcmin^ of the GOODS-S field was observed with 60s frame times. Most of the area 
within the GOODS fields had four dithers per position while the UDF area had 32 dithers 
per position. Each frame was distortion corrected and background subtracted. The dithered 
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frames were combined using "drizzle" (IFruchter fc Hookll2002l ). Due to the limited number 



of dithers per position, point kernel drizzling, which would have minimized the correlated 
noise in the mosaics was not possible. The final mosaics have a plate scale of 0.9" per pixel. 
The exposure time in the final mosaics is ~4 min pix~^ over most of the GOODS area while 
the UDF area was observed with 32 min pix~^. JD2, although in the UDF region, seredip- 
itously lies in a part of the coverage map where the exposure time is 68 min, twice that of 
the nominal UDF. As a result, the statistical uncertainty on its fiux is lower than for the 
average UDF source of the same brightness. 

The spatial resolution of Spitzer at 16/im and 22/im is 4.1" and 5.2" full width at half 
maximum (FWHM), respectively. At this resolution, the majority of sources, including JD2, 
are point sources. This makes it possible to use prior positions from higher resolution data, 
such as the GOODS 3.6— 8 /im imaging (jDickinson et al. 2003 . Dickinson et al., in prep.) 



and apply point-source fitting techniques to measure the fiux density of the source. The 
technique has successfully been used for cataloging 24/im sources from the GOODS/MIPS 
imaging survey. The results are consistent with aperture photometry, with appropriate 
aperture corrections, for isolated sources and alleviate the contamination to the photometry 
from the wings of the point spread function for sources which have nearby companions. 
Postage stamp image cutouts of JD2 are shown in Figure 1. 

The 16/im and 22/xm fiux densities of the source are 13.5±3.5/iJy and 56±12/iJy re- 
spectively. In comparison, the 24 fim fiux density of the source from the GOODS imaging of 
the field is 51±4yuJy. HUDF-JD2 is separated by 7" from a brighter mid-infrared source to 
the south east which is a z = 0.457 spiral galaxy (Stern et al., in prep.). In order to estimate 
the systematic uncertainty associated with fitting the fiux density of a faint source near a 
bright one, we performed a Monte-Carlo simulation. We identified a relatively isolated source 
in the mosaic with about the same fiux density as the spiral galaxy: this is the galaxy at 
3:32:43.49, -27:45:56.45 (J2000) which has a fiux density of 121±10/iJy at 16/im. We added 
an artificial source at a distance between 7" and 8" from this source and extracted the fiux 
density using the positional priors. The process was repeated 100 times each for a range 
of fiux densities straddling the measured fiux density of the source. The extracted fiux was 
compared with the input fiux density to assess the systematic uncertainty and/or fiux bias 
in the measurement. For a source with a brightness comparable to JD2, the Monte-Carlo 
analysis yielded la uncertainties of 14/iJy, 12/iJy and 10/iJy at 16, 22 and 24 /im, respec- 
tively. This systematic uncertainty is primarily due to the extended wings of the brighter 
source and is a factor of ~2-3 larger than the statistical uncertainty. We therefore adopt 
values of 14±14/iJy at 16/im, 56±12/iJy at 22/im and 51±10/iJy at 24/im for JD2. 



The sharp increase in fiux density between 16/im and 22/im can either be due to the 
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entrance of polycyclic aromatic hydrocarbons into the 22 and 24 /xm bandpasses or the 9.7/im 
sihcate absorption feature entering the 16/im bandpass. A third possibihty is the decreasing 
contribution from hot dust surrounding an AGN because the dust is close to its subhmation 
temperature. We evaluate each of these possibilities in the following sections. 

We note that the IRS is calibrated with respect to a constant z/F,y source spectrum. Color 
correction terms to the observed flux are smaller than 1% for the 16/im and 22/im data for 
the observed source spectrum and can be neglected. Furthermore, when we estimate fits to 
the observed photometry, we integrate template spectra over the filter bandpasses (Figure 
2). 



3. Derived Source Parameters 



The o ptical/ near-inf r ared photometry of JD2 has been discussed in iMobasher et al. 



(120051 ) and lDunlop et al.l (120071 ). Spectral energy distribution fits to the photometry, which 
provide a photometric redshift solution, show a bimodal distribution of minimum values 
at z ~ 2.5 and z ^ 6.5. On e set of value s are f avored over the other depending on the 
choice of optical photometry. iMobasher et al. J2005h. with th e ir stri ngent BViz limits, derive 
the probability of a z < 5 source to be 15%. iDunlop et al.l (120071 ) instead adopt marginal 
detections of the source in the Viz bandpasses which we discuss in Secion 3.2. Any significant 
detection of the gal axy at optical wavele ngths would rule out a photometric redshift solution 
which favors z > Q. iDunlop et al.l (120071 ) find a minimum in the distribution at 2; ~ 2.2. 



The detection of HUpF-JD 2 in the deep GOODS 24/im images did not break the de- 



generacy. IMobasher et al.l (120051 ) suggests the hot dust emission is from an obscured AGN 
in the galaxy. Alternately, the 24 /xm emission could be dominated from polycyclic aromatic 
hydrocarbon (PAH) feature s which would be present if the ob ject is an extremely red, dusty 
starburst galaxy at z ~ 2.5 (lYan et al.ll2004l : IStern et al.ll2006l ). Inclusion of the 16 and 22 /xm 
data allows us to measure a photometric redshift from the mid-infrared which could poten- 
tially break this degeneracy, independent of the optical/near- infrared SED fits. Although 
fitting optical to mid-infrared simultaneously is more elegant, it is virtually impossible since 
there is virtually no correlation between optical/near-infrared SEDs and mid- infrared SEDs 
of galaxies. 
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3.1. Mid-Infrared Photometric Redshifts 

The mid-infrared SED of a star-forming galaxy is a complicated interplay of warm dust 
continuum, PAH emission and silicate absorption features. The ratio of flux densities in 
the three abutting 5^)«teer passbands at 16, 22 and 24 /xm allow redshifts of objects to be 
constrained as these dust features move through these filters. PAH are only 0.5-1/im wide, 
while the 9.7/xm silicate feature is 3— 4/im wide, depending on the continuum level adopted. 
As a result, adjacent passbands can show widely different flux density ratios as a function 
of redshift (Figure 3). 

At 0.2 < z < 0.6, the 9.7/im silicate feature enters the 16/im bandpass while the 22 
and 24 /im bandpasses are tracing the warm dust and 11 and 12 /im PAH features. At these 
redshifts, sources show low 16yum/22/im flux density ratios. At 0.6 < z < 1.2, the 9.7/im 
features moves out of the 16 /xm bandpass and into the 22 /im bandpass, while the 6.2 /im 
and 7.7/im PAH features fall in the 16 /im window, boosting the 16/im/22/im ratio. At 
1.2 < z < 1.7, the 9.7/im feature moves into the 24 /im band. Since the 22 /im bandpass is 
~2 /im wider than the MIPS 24 /im band, the 7.7/im feature enters the 22 /im bandpass. The 
net effect is to boost the 22/im/24/im flux ratio while lowering the 16/im/22/im ratio. At 
z > 1.9, the PAH emission shifts out of the 16 /im band, decreasing the 16/im/22/im ratio 
while the 22/im/24/im flux ratio varies as the 6.2 and 7.7/im PAH features move within 
the passbands. This variation in flux ratios is illustrated in Figure 3 for four different SEP 
types derived frona spectral observations of galaxies in the local Universe (lArmus et al.ll2007l : 
Brandl et aDl2004h . 



While it is true that variation of PAH line ratios or strength of silicate absorption can 
cause broadband fluxes to vary signiflcantly, the combination of photometric redshifts from 
the optical/near-infrared and mid- infrared photometric redshifts using these three band- 
passes can help break redshift degeneracies. 

The flux density ratios observed for HUDF-JD2 are a 16/im/22/im ratio of 0.25" 



-0.19 
-0.08 

and a 22 /im/24 /im ratio of l.llg'^g. This indicates that the galaxy must either be at 2; ~ 0.6 
if it has very strong silicate absorption or at 2; ~ 1.7 if it is a typical starburst galaxy with 
strong PAH emission. 

In Figure 4 and Table 1, we illustrate the quality of flts to the photometry at 16, 
22 and 24/im from different mid-infrared template spectral energy distributions. We con- 



PAH and strong silicate absorption ( 


Charv & Elbaz 




2001 




. ■i'orster-Schrei 


3er et al. 


2001; 


Le Floc'h et al. 


2001; 


Dale & Helou 


2002: 


Brandl et al. 


2004; 


Armus et al. 


2007 


). We also 



consider the non-detection of JD2 in deep 70 /im observations of the GOODS flelds which 
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achieve 5a flux density limits of ~2 mJy (IFrayer et al.ll2006l ). The best fit to the 16/im, 22 /im 
and 24 fim photometry is from a starburst source with strong PAH emission at 2; = 1.7. The 
increase in the 22/im and 24/xm fiux relative to 16yum is due to the entrance of the 7.7/im 
PAH complex into the 22/im and 24/im bandpasses. The fiux in the 16/im bandpass is lower 
since the 6.2/im PAH feature is much weaker. 

Extrapolating the templates that are fit to the mid-infrared photometry also yields an 
estimate of the far- infrared luminosity for the galaxy. The Lir=L(8— 1000 fim) luminosity of 
the source is ~5xl0^^ Lq, which corresponds to a star-formation rate of 80 Mq yr~^. We 
compare this with optica l/UV estirnates o f star-formation in the following section. Using the 
radio-FIR correlation in lYun et al.l (l200ll ). we predict the 1.4 GHz to be 12/iJy and the 8.4 
GHz fiux density to be 3.7/iJy while the predicted fiuxes at 70 fim and 850 fim are shown in 
Table 1. The deepest 70 fim and 850 /xm surveys are currently sensitive to 1.5—2 mJy. The 
predicted intensities of JD2 are below these limits implying that confirmation of the derived 
far-infrared luminosity will have to await deeper observations by Herschel and ALMA. 

values are also reduced for an obscured AGN template aX z = 0.6 which has the 
observed 16 fim fiux density suppressed relative to the fiux density at 22 fim and 24 fim due 
to the entrance of the 9.7 fim silicate absorption feature. An obscured AGN template like 
Mrk231 can also fit the 8/im fiux density of the source as shown in Figure 4 while the 
starburst template fits at z ~ 1.7 require the optical/near- infrared fiux to be dominated 
by starlight (Figure 6). The no n-detection of any spe ctral lines in the optical/near-infrared 
spectroscopic data presented by lMobasher et al.l (120051 ) . as well as the difficulty in accounting 
for the non-detection of the source at optical wavelengths (Section 3.3) makes the z ~ 0.6 
hypothesis unlikely. On the other hand, z ~ 1.7 lies within the so-called redshift "desert", 
where the [OH] A3727 doublet, [OIIl] A5007 line and Hq;+[NII] lines are all severely affected 
by atmospheric transmission and strong OH sky lines. Furthermore, as we discuss in Section 
3.3, at z = 1.7, the red SED of the source between 2.2/im and IRAC 3— 8/im could be due 
to the broad 1.6 fim bump in the stellar SED being in the IRAC passbands. 



m 



We now consi d er the implications of the IRS detections to the high redshift advocated 
Mobasher et al.l (120051 ). If we adopt the hypothesis that JD2 is at z ~ 6.5 and the 24 fim 



emission is dominated by AGN light, we must first subtract t he stellar contrib u tion t o the 
mid-infrared fiux densities. We find that for the SED fit by iMobasher et al.l (120051 ). the 
contribution of stellar photospheric emission to the fiux densities at 16, 22 and 24 fim are 
10.7, 6.5 and 6.0 fiJj respectively. We then attempt to fit the starlight- subtracted mid 



i nfrar ed photometry with the two mid-infrared AGN templates presented in iMobasher et al. 
((20051). The values are significantly worse than the best fits at z ~ 1.7 (Table 1 and Figure 



5). The Mrk231 template, fairly typical of an obscured AGN, is bluer in its 16 /im/22 fim fiux 
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ratios at z = 6.5 than what the observations indicate. NGC1068, which is a Compton-thick 
AGN with decreasing hot dust continuum at shorter wavelengths due to dust subhmation, 
would best fit the observed mid- infrared photometry at z = 3.7. The hot dust emission falls 
off too rapidly for it to account for the 22 yum flux density with respect to the 24 fim flux 
density if it were at z = 6.5. One could imagine that varying the AGN template could allow 
for the mid- infrared photometry here to be fit but given that the 2; = 1.7 fit using typical 
mid-infrared templates agrees with the parameters derived from the optical/near-infrared 
fits, we find the z = 6.5 interpretation for JD2 substantially weakened. 



3.2. The Optical/Near-infrared Photometry 

Given the derived redshift based on the mid-infrared photometry of the source, we 
critically investigate the optical/near- infrared photometry to assess potential sources of error. 

The ACS BViz photometric measurements are very important to the interpretation 
of JD2, since a significan t optical detection wo uld almost certainly rule out the z ~ 6.5 
hypothesis. We note that iMobasher et al.l (120051 ) erroneously reported that their ACS pho- 
tometry was measured in an 0'.'9 diameter aperture, when in fact a smaller, 0'.'48 diameter 
aperture was used (M. Dickinson, private communication). A larger aperture size would lead 
to more conservative photometric limits, and is probably to be recommended, given the size 
of the galaxy as measured from the NICMOS if-band images (measured half-light radius 
0'.'3, uncorrected for PSF effects, see Mobasher et al.). However, a larger aperture risks 
including light from three faint, neighboring galaxies located about 1" away from JD2. We 
have masked out regions with diam;t« a'9 around tho=e three gala^iefl before measuring 
photometry for JD2 in a I'.'O diameter circular aperture. 



The ACS images were drizzled using a point kernel (see iBeckwith et al.l |2006| ) which 
should lead to noise that is uncorrelated between adjacent pixels. We verified that this is 
the case by measuring the autocorrelation function of the noise after maski ng galaxies. We 



measured the noise on various scales, including 50 pixel apertures as used by IBeckwith et al. 



(120061 ). Our mos t conservative ( l argest ) noise measurements are 31 to 37% smaller than 
those reported by IBeckwith et al.l (l2006l . their Table 5). 



Dunlop et al.l (120071 ) report faint detections of positive fiux in the ACS V, i and z- 



bands, which, taken together, drive their photometric redshift estimate to favor a lower 



^In another typographical error. iMobasher et ahl (|2005l ) reported masking regions 0'.'5 in diameter around 
the neighboring galaxies. In fact, the masking diameter was 0'.'9, which we also adopt here. 
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value, z ^ 2.15. We do not reproduce these measurements; using a I'.'O diameter aperture 
and our noise measurem ents, and the c o rrecti on of ~0.1 mag for the energy falhng outside 
the aperture reported in ISirianni et aLl (120051 . Table 3), we find the 1o photometric limits 
given in Table 2. These are ~0.8 mag brighter than the limits reported in Mobasher et 
al. (2005), mainly due to the larger aperture adopted here. It is possible that the positive 
flux detected in the Dunlop et al. measurements arises in part from the faint neighboring 
galaxies. 



The NICMOS Jno photometry and ifiBn photometry i n iMobasher et al.l (l2005l ) are very 
similar to those in the UDF catalog of [Thompson et al. J2005h . As a res ult, we have not 



remeasured these values. However, we note that the photometric values in iMobasher et al. 



m 



(I20051) do not in c lude t he NICMOS count rate dependent non-linearity correction discussed 



de Jong et al.l (120061 ). The effect of this non-linearity correction is to make the photometry 
brighter by 0.22 mag in the Jno band and 0.12 mag in the -ffieo band as shown in Table 2. 

We also remeasured the photometry of the source at IRAC wavelengths. We used SEx- 
tracto r catalogs and measured the photometry in 3" diameter beams, similar to IMobasher et al, 
( I2OO5I ). To derive aperture corrections and systematic uncertainties, we adopt a different ap- 
proach. The primary systematic uncertainty in this measurement is the uncertain sky level 
due to the bright galaxy 7" to the South-East. We input an artificial point source whose 
brightness is that of JD2 into the final mosaic. The centroid of the source with respect to 
the bright galaxy is kept similar to that of JD2, i.e the artificial source has a centroid which 
is between 7" and 8" from the galaxy. The SExtractor routine was run on this fake image 
and the photometry of the fake source measured. We find a systematic uncertainty whereby 
the extracted flux of the fake source in the catalog was brighter at 3.6 and 4.5 \xm but fainter 
at 5.8 and 8.0 /xm. The process was repeated a 100 times to measure the average systematic 
and statistical error using the corresponding point spread function at each of the four IRAC 
wavelengths. We revise the aperture corrected magnitudes for JD2 to 22.24, 22.00, 21.63 
and 21.51 AB mag at 3.6, 4.5, 5.8 and 8.0 /im respectively. There is, however, a correction 
to be applied to these magnitudes. 

There are three faint galaxies within 1" of JD2. Two of these are rather faint and blue 
while the brightest of the three falls within the IRAC 3" beam. These galaxies are detected at 
all wavelength between the B and ifieo band. To estimate the contribution of these sources 



to the photometry within the IRAC beam, we fit iBruzual &: Chariot I (120031 ) templates to the 
multiband photometry of these sources, leaving redshift as a free parameter. We use the best 
fit SEDs to estimate the redshift and corresponding flux densities at the IRAC wavelengths. 
We find the brightest of these three sources, located at 3:32:38.76,-27:48:28.91 (J2000), to be 
at redshift 3.4 with magnitudes of 25.95, 25.83, 25.82 and 25.80 at 3.6, 4.5, 5.8 and 8.0 /xm 
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respectively. The second source, located at 3:32:38.81,-27:48:39.79, is at z ~ 3.1 with a 
contribution that is about 0.6 mag fainter at each of the four passbands. We subtracted the 
contribution of both these sources from the photometry of JD2. 

To summarize, the corrected limits for JD2 in the ACS bands are about 0.8 mag brighter 
due to the larger aperture, the photometry in the NICMOS Jno and Hiqq bands is brighter 
by ~ 10 — 20% due to the non-linearity correction and t he photometry in the IR AC bands 
is about 10% fainter on average than those adopted in iMobasher et al.l (120051 ). The net 
effect is to reduce the amplitude of the apparent near-infrared to IRAC "break" slightly. 
The optical/near-infrared photometry of JD2 after these corrections have been applied is 
shown in Table 2 and which we fit for in S ection 3.3. Fo r the s ake o f completeness, w'e also 
investigate redshift constrained fits to the iDunlop et al.l (120071 ) and iMobasher et al.l (120051 ) 
optical/near-infrared photometry. 

In addition, we attempted to fit elliptical isophotes to the galaxy to the SE and sub- 
tracted it out from the image. This process left significant residuals due to the asymmetric 
nature of the Spitzer/IRAC point spread function. The aperture corrected point source pho- 
tometry we measure for JD2 in images with the galaxy subtracted is 22.36, 21.95, 21.72 and 
21.64 AB mag. The difference between these values and those quoted in Table 2 can be 
attributed almost entirely to the difficulty in measuring the absolute sky background. Thus, 
we provide photometric uncertainties of 0.15 mag for all 4 bands. 



3.3. Revisiting the Optical/Near-Infrared Photometric Fits 



We first evalu ate the quality of the fits to the optical/near-infrared photometry at 
z = 0.6 using the iBruzual fc Chariot I (120031 ) population synthesis models. The range of 
parameter space for the models was solar and 0.2 solar metallicity, 19 e-folding timescales 
(r) for the starburst from an instantaneous starburst to constant star-formation, ages (t) 
from 0.1 Myr to the age of the Universe at the chosen redshift and extinction (Ay) between 
and 10 mags. At ^ = 0.6, the 8 /im fiux density is affected by dust emission and is ignored 
in the fits. 

If JD2 is a 2 = 0.6 starburst like NGC6240 but with silicate absorption which sup- 
presses the observed 16 /xm fiux density, the weakness/non-detection of the source in the 
BV iz passbands require s A y=9.9 mag of visual ex tinction irrespective of whether we adopt 
the iDunlop et al.l (120071 ) or IMobasher et al.l (120051 ) photometry. This is not impossible but 
unusually large. It would imply a young 8 Myr old starburst with a mass of 2x10^ M© and 
a true star-formation rate, as derived from extinction correction of the ultraviolet fiux, of 
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~100 M0 yr~^, all of which is dust obscured. The fits to the mid-infrared photometry at 
z = 0.6 imply a much lower dust-obscured star-formation rate of ~5 Mq yr~^. Although it 
is possible for optically thick star-formation, such as that found in local ULIRGs, to result 
in the mid- and far-infrared luminosity of galaxies exceeding their extinction corrected ul- 
traviolet luminosity, the converse is rarely true. This weakens the possibility of a z = 0.6 
star burst. 

If the source is a 2; = 0.6 obscured AGN like Mrk231, then the contribution of the hot 
dust around the AGN to the near-infrared photometry needs to be subtracted before fitting. 
For the Mrk231 template, we find this to be 0.8, 1.2, 1.4, 2.6, 3.5, 5.6 and 10.5//Jy at 1.2, 1.6, 
2.2, 3.6, 4.5, 5.8 and 8.0/im respectively. The photometry at 1.1 and 1.6//m of the source 
is a factor of 3 — 10 below the expected values again indicating that a Mrk231 type source, 
while consistent with the mid-infrared photometry, cannot fit the entire SED. 

Thus, a z = 1.7 dusty star- forming galaxy SED is most consistent with all the pho- 
tometry. We adopt this redshift and re-fit the optical/near-infrared data using the BC03 
models. Results are shown in Table 3. HUDF-JD2 is best fit with a solar metallicity tera- 



plate from the BC03 library with Ay ~ 4 mag of visual extinction using a ICalzettil (120011 ) 
dust extinction law. We caution that the metallicity is not strongly constrained. The best-fit 
model is a post-starburst BC03 SED with t = 570 Myr, r= 30 Myr and a st e llar m ass of 



6x10^° Mq. This is similar to the 2; ~ 2.5 fits presented in iMobasher et al. J2005h . The 



result is not surprising since post-starburst galaxies could have significant dust content from 
dust produced in supernovae and AGE stars. In fact, it has been observationally shown 
that the host galaxies of core-collapse and Type la supernova e are dustier than field galaxies 



with the same observed optical brightness (IChary et al.ll2005l ). We caution however, that the 



post-starburst hypothesis, motivated by the large value of t/r, would be weakened if most 
of the star-formation takes place in optically thick regions of the galaxy. 

From the SED fits, we estimate that the star-formation rate, as measured by the 1500A 
flux escaping the galaxy, is very small (~2xl0~^ M© yr~^). Application of an extinction 
correction to the UV luminosity results in a star-formation rate of <1 M© yr~^. In contrast, 
its true star- formation rate as derived from the mid- infrared fits is ~80-90 M© yr~^. The 
fact that JD2 is a compact near-infrared source with a scale length of ~0.3", combined with 
the large ratio of infrared to ultraviolet star-formation rate estimates, strongly suggests that 
JD2 is being powered by a compact, optically-thick starburst in its nucleus. 



There is a caveat. If we use the iBruzual fc Chariot I (120031 ) models to solve for an 



independent redshift from the optical/near-infrared photometry presented in this paper, 
we continue to get a high redshift minimum at z = 7.2 which has a reduced of 1-05 
compared to the reduced values of 1.9 and 2.4 ai z = 2.4 and z = 1.7 respectively. Thus, 
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optical/near-infrared photometric redshifts continue to favor the high redshift hypothesis, 
while mid-infrared photometric redshifts favor the lower redshift hypothesis. We favor the 
mid-infrared photometric redshifts derived here since optical photometric redshifts for dusty 
galaxies have been shown to have a large scatter based on a comparisoi i between spe ctroscopic 
and photometric redshifts for 24 /im detected sources in GOODS-N (jCharyll2006l ). 

A comparison of the prop erties of HUDF- j p2 w' ith the massive, red galaxy population 
at 1.5 < 2; < 3 is illustrative. iPapovich et al.l (120061 ) show that the star-formation rate of 
galaxies scale with the stellar mass at these redshifts and that the specific star-formation 
rate of a dusty starburst with mass of ~6xl0^° Mq is about ~ 1 — 20 Gyr~^. UDF-JD2, 
for z ~ 1.7, has a specific star-formation rate of 1.6 Gyr^^ which is consistent with the 
massive, dusty starbursts at these redshifts. The optical extinction derived for this galaxy 
is rather large, corresponding to E{B — l^)=1.2. Spectroscopically confirmed galaxies at 
these redshifts which are 24/im detected have E{B values derived from their ultraviolet 
slopes of ~0.3. This is most likely a bias since spectroscopic confirmation requires the 
detection of optical/ultraviolet lines which are preferentially detected in galaxies that are 
relatively transparent and have low extinction values. A bet ter comparisoii is wi th the 
derived extinction values for local dusty starbursts presented in iHopkins et al.l (I2OOII ). The 
average E{B — to the stellar continuum for that sample at an infrared luminosity of 
5x10^^ Lq is about 0.25. Thus, HUDF-JD2 does not appear to be dissimilar from dusty, 
luminous infrared galaxies (LIRGs) at 2; ~ 2. However, its derived E{B — appears to 
be a factor of ~4 higher than other LIRG samples, consistent with its significant detection 
at mid- infrared wavelengths and compact morphology in the NICMOS 1.6 yum image which 
suggests a compact, nuclear starburst similar to many local infrared luminous galaxies. 

JD2 is about 1.4 mags redder in observed i?— [3.6] colors and about 3 r nags fainter 
at 3. 6 /im than HRIO, a 2; = 1.44 spectroscopically confirmed dusty starburst (jStern et al. 



20061 ). The predicted Ha line flux for this galaxy, including the extinction in the line, is 
~10~^^ ergs~^ cm~^. This is a factor of 3 below the sensitivity of the existing observa- 



tions in 



spectrograph. 



Mobasher et al.l (120051 ) but should be detectable with a space-based near-infrared 



4. Conclusions 

We present new 16 and 22 /xm mid-infrared photometry of the HUDF source JD2 which 
had tentatively been classified as a z ~ 6.5 massive galaxy. The mid-infrared spectral energy 
distribution shows a sharp increase in flux density between 16 /xm and 22 /xm, most consistent 
with the 7.7/im PAH feature entering the 22 /im bandpass. Our analysis concludes that the 
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galaxy is a starbursting, luminous infrared galaxy (LIRG) at a redsliift of z ~ 1.7. The 
non-detection of the source in the UDF BViz data is due to the presence of Ay ~ 4 mag of 
dust internal to the galaxy. We have re-e yaluated the optical /ne ar- infrared photometry of 
this source and find the stringent limits in iMobasher et al. I J2005h are more consistent with 
the data rather than the marginal detection reported in iDunlop et al.l (120071 ). Photometric 
redshift fits to the optical/near-infrared data alone continue to favor a z > 6 solution. 
However, the unusually high Lmir/Lx for JD2 compared to other obscured AGN, the large 
scatter in optical photometric redshift solutions for dusty galaxies and the similarity between 
the multiwavelength properties of this galaxy and other dusty starbursts in the 1.5 < z < 3 
range lead us to favor the mid-infrared photometric redshift solution. The non-detection of a 
z 2 LIRG in the deepest optical data taken to date strongly cautions against interpreting 
near-infrared dropout sources in shallower surveys as 2; > 6 galaxies. 
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Fig. 1. — 16, 22 and 24/im snapshots of the Hubble Ultradeep Field source JD2 (circled) 
illustrating the quality of Spitzer mid-infrared data on this source. Images are ~40" on a 
side with North up and East to the left. JD2 is at 3:32:38.74, -27:48:39.9 (J2000). 




Fig. 2. — Response curves for the IRS 16 /xm and 22 /xm bandpasses and MIPS 24 /xm band- 
pass. 
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Fig. 3. — 16;um/22 /iin and 22 /im/24 yum flux density ratios as a function of redshift for 
typical infrared luminous galaxies at 0< z < 3. The sources shown are NGC7714 - a typical 
starburst (SB) with strong PAH, NGC6240 - a composite AGN+SB object with strong PAH, 
Arp220 - a starburst with strong silicate absorption and weak PAH and Mrk231 - an AGN 
with weak silicate absorption. The range of flux ratios of HUDF-JD2 with statistical errors 
is shown by the horizontal dot-dash lines. The flux ratios indicate a starburst with strong 
PAH at 2; ~ 1.7 or a source with strong silicate absorption at z ~ 0.6. 



Table 1. Results of Mid-infrared Spectral Energy Distribution Fits 



Template Adopted 


Best fit redshift 






Predicted Fluxes 


Type of Object^^ 


Notes 






L0 




(70, 850Atm in mjy) 






Dale &; Helou 


1.7 


4.7x10" 


0.30 


0.6, 0.5 






Chary & Elbaz 


1.6 


5.6x10" 


0.37 


0.9, 0.4 






NGC7714 


1.6 


2.5x10" 


0.42 


0.4, 0.2 


Starburst 




M82 


1.7 


1.6x10" 


0.48 


0.2, 0.09 


Starburst 




NGC6240 


0.5 




0.60 


1.0, 0.02 


Compton thick AGN+SB 


large Ay 














inconsistent UV/IR SFR 


Mrk231 


0.5 


6.9x10^5 


0.76 


0.4, 0.004 


Obscured AGN 


exceed NIR fiux 


Arp220 


1.8 


2.1x10^2 


0.95 


1.7, 1.5 


Starburst 




NGC1068-nucleus 


3.7 


1.4x10^2 


1.02 


0.3, 1.0 


Obscured AGN 


Worse 


Mrk231 


6.5 


7.4x10^2 


2.0 


0.2, 2.4 


Obscured AGN 


Worse 


NGC1068-nucleus 


6.5 


1.3xl0i-^ 


1.1 


0.6, 7.1 


Obscured AGN 


Worse x^ 



''For the star-forming templates, the star-formation rate is 1.71x10 ""^"xLir. 

''For the AGN templates at 2 = 6.5, the stellar photospheric emission has been subtracted from the mid-infrared photometry. To identify the 
best-fit mid-infrared redshift for the AGN templates the actual mid-infrared photometry was used. 
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Fig. 4. — (Left panel) Spectral energy distribution fits to the mid-infrared photometry 
of HUDF-JD2 along with the distribution of values with redshif t (Right panel) . The 
SEDs are plotted at their best fitting redshift which is z = 1.7 for the lDale fc Heloul (120021 ) 
starburst templates, z = 0.5 for NGC6240 and z = 0.5 for Mrk231. The preferred fit (i.e. 
lowest value) indicates the source is a luminous infrared galaxy with strong PAH emission 
at 2 ~ 1.7. A low is also obtained at z ~ 0.6 from a source with either strong silicate 
absorption or the gap between the 8 and 11 /im PAH complexes. However, the 2; ~ 1.7 fit is 
more consiste nt with the fits to th e optical/near-infrared photometry. The stars in the left 
panel are the iDale fc Heloul (120021 ) SED convolved with the filter curves. 
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Observed X [Atnn] Redshift 



Fig. 5. — (Left panel) Fits to the mid- infrared photom e try of HUDF-JD2 at z = 6.5 using 
the two AGN templates presented in iMobasher et al.l (l2005l ). The solid red and purple 
circles are the AGN templates of the corresponding color integrated through the relevant 
bandpasses. The right panel shows the distribution of values. The solid black square 
shows the minimum obtained at z = 1.7 for the starburst spectral energy distribution 
discussed in Figure 4 and in Section 3.1. If the source were really at z = 6.5, there would be a 
stellar contribution to the mid- infrared photometry as discussed in the text. Although, this 
improves the formal estimates, the z = 6.5 fits still yield substantially worse estimates 
than the best fit at 2; = 1.7 as shown in Table 1. 
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Fig. 6. — Optical to mid-infrared spectral energy distribution of HUDF-JD2. Solid black 
square s are the photo n aetry presented in this paper. Solid red squares are the photome- 



trv in iMobasher et al.l (I20051) . The optical/near-infrared photometry has be en fit with the 



Bruzual fc Chariot I (120031 ) model and the mid- infrared photometry fit with the lDale fc Helou 



(120021 ) model. The source is most likely a luminous infrared galaxy (LIRG) with strong PAH 
emission at 2; ~ 1.7. The inferred star-formation rate is ~80 Mq jt~^ with ~4 mag of visual 
extinction. 
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Table 2. Revised Optical/Near-infrared Photometry for JD2 



Wavelength Photometry (AB mag)'' 

B 
V 
i 
z 

Jiw 
HiQo 
Ks 

3.6 /xm 
4.5 /xm 
5.8 /im 
8.0 /im 

16 lira 
22 /jm 
24 /im 



^Photometric Hmits, where provided, are 
2(T. See text for details. 



Table 3. z=1.7 SED Fits to Optical-Near-infrared Photometry 



Mass 


Av 


Age 


r 


UV SFR 


Corrected SFR 




Notes 


M0 


mag 


Myr 


Myr 


Mq yr-i 


Mq yr-i 






6.3xl0io 


3.6 


570 


30 


1.7x10"'' 


0.74 


2.4 


Photometry from Table 2, Solar metallicity BC03 


9.8x1010 


4.7 


400 





8x10-5 


5 


3.9 


Mobasher et al. photometry, Solar metallicity BC03 


1.3x10" 


4.6 


1020 


30 


4x10-5 


1.8 


4.0 


Mobasher et al. photometry, 0.2 solar metallicity BC03 


4.6x10^° 


4.3 


453 


700 


0.003 


93 


1.6 


Dunlop et al. photometry. Solar metallicity BC03 


1.9xl0i° 


4.7 


47.5 


50 


0.005 


308 


1.6 


Dunlop et al. photometry, 0.2 solar metallicity BC03 



^^4 degrees of free dom with 7 data points and 3 parameters (Age, Extinction and r) for the iMobasher et al] l|2005l ) photometry. The 
iDunlop et al. I ll2007t) photometry has 10 data points. 



>29.7 

>30.1 

>30.0 

>29.3 

26.8±0.3 

24.82±0.07 

23.95±0.132 

22.29±0.15 

22.05±0.15 

21.67±0.15 

21.54±0.15 

21 n4+^'*® 

19.53±0.21 
19.63±0.21 



